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Abstract 

Q_i' Possible models of Yukawa interaction are discussed in the two Higgs doublet model (THDM) 

under the discrete symmetry imposed to avoid the flavor changing neutral current at the leading 
. order. It is known that there are four types of such models corresponding to the possible different 

> 

I/*} ' assignment of charges for the discrete symmetry on quarks and leptons. We first examine decay 

' 

SO ■ properties of Higgs bosons in each type model, and summarize constraints on the models from 

. current experimental data. We then shed light on the differences among these models in collider 

o. 

Q^ phenomenology. In particular, we mainly discuss the so-called type-II THDM and type-X THDM. 

o 

The type-II THDM corresponds to the model with the same Yukawa interaction as the minimal 
^ ■ supersymmetric standard model. On the other hand, in the type-X THDM, additional Higgs 

h : 

bosons can predominantly decay into leptons. This scenario may be interesting because of the 
motivation for a light charged Higgs boson scenario such as in the TeV scale model of neutrino, 
dark matter and baryogenesis. We study how we can distinguish the type-X THDM from the 
minimal supersymmetric standard model at the Large Hadron Collider and the International Linear 
Collider. 

PACS numbers: 12.60.Fr, 14.80. Cp, 14.60.-z 
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I. INTRODUCTION 



The CERN LHC, which is a proton-proton collider with maximal energy 14 TeV, has 
started its operation 1|. The most important purpose of the LHC is hunting for the Higgs 
boson, the last unknown particle in the standard model (SM). The Higgs sector is introduced 
for the spontaneous breakdown of the electroweak gauge symmetry. Weak gauge bosons (W 
and Z) then receive their masses via the Higgs mechanism, and quarks and charged leptons 
receive their masses through the Yukawa interaction. Because the SM Higgs sector has not 
been confirmed yet, the possibility of its non- minimal form should also be considered in 
order to understand the nature of the symmetry breaking. In fact, it is well known that 
many candidates of physics beyond the SM predict extended Higgs sectors. For example, 
super symmetric extensions of the SM have at least two Higgs doublets [2] . Successful models 



based on dynamical symmetry breaking a 
sector in the low energy effective theory 



so may require the non-minimal form as the Higgs 
^j]. In addition, some of new physics models that 



are intended to solve the problems which the SM cannot explain, such as tiny neutrino 
masses, essence of dark matter, and baryogenesis, has been built with the extension of the 
electroweak symmetry breaking sector (4, 5, 6, 7|. 

There are two basic experimental constraints which an extended Higgs sector must re- 
spect: those on the electroweak rho parameter (p) as well as the flavour changing neutral 
current (FCNC). The measured value of the rho parameter (p exp ~ 1) suggest that the 
electroweak sector of the model would approximately have a global SU(2) symmetry (the 
custodial symmetry) which guarantees p = 1 at tree level 8j. For this to be satisfied, a 
possibility of the structure with (multi-) isospin doublets (plus singlets) would be natural 
as the Higgs sector. In the SM, FCNC phenomena are suppressed due to the electromag- 
netic gauge symmetry and due to the Glashow-Iliopoulos-Maiani mechanism [9] , so that the 
experimental bounds on the FCNC processes are satisfied. In models with more than one 
Higgs doublet, this is not true in general, because two or more Yukawa matrices for each 
fermion cannot be simultaneously diagonalized. It is well known that, to avoid such Higgs- 
boson-associated FCNC interactions, each fermion should couple to only one of the Higgs 
doublets. This can be realized in a natural way by imposing a discrete Z2 symmetry 10]. 

In this paper, we discuss phenomenological differences in various models of Yukawa in- 
teractions under the discrete Z 2 symmetry in the two Higgs doublet model (THDM). It has 



2 



been known that there are four patterns of the Yukawa interaction depending on the assign 



0Q- 



We refer them 



ment of charges for quarks and leptons under the Zi symmetry [U, 
as type-I, type-II, type-X and type-Y THDMs 1 . The type in the Yukawa interaction can be 
related to the new physics scenarios. For example, the Higgs sector of the minimal super- 
symmetric standard model (MSSM) is the THDM with a supersymmetric relation [2j among 
the parameters of the Higgs sector, whose Yukawa interaction is of type-II, in which only 
a Higgs doublet couples to up-type quarks and the other couples to down-type quarks and 
charged leptons. On the other hand, a TeV scale model to try to explain neutrino masses, 
dark matter and baryogenesis has been proposed in Ref. 7J. In this model the Higgs sector 
is the two Higgs doublets with extra scalar singlets, and the Yukawa interaction corresponds 
to the type-X, in which only a Higgs doublet couples to quarks, and the other couples to 
leptons. Therefore, in order to select the true model from various new physics candidates 
that predict THDMs (and their variations with singlets), it is important to experimentally 
determine the type of the Yukawa interaction. 

There have been many studies for the phenomenological properties of the type-II THDM, 
often in the context of the MSSM [2]. On the contrary, there has been fewer studies for the 
other types of Yukawa interactions in the THDM. The purpose of this paper is to clarify 
phenomenological differences among these types of Yukawa interactions in the THDM at the 
LHC and the International Linear Collider (ILC) [3]. We first study the decay rates and 
the decay branching ratios of the CP-even (h and H) and CP-odd (A) neutral Higgs bosons 
and the charged Higgs bosons {H ) in various types of Yukawa interactions. It is confirmed 
that there are large differences in the Higgs boson decays among these types of Yukawa 
interactions in the THDM. In particular, in the case where the CP-even Higgs boson h is 
approximately SM-like, H and A decay mainly into t + t~ in the type-X scenario for the wide 
range of parameter space, while they decay mainly into bb in the type-II scenario. We then 
summarize constraints on the mass of H from current experimental bounds in various types 
of Yukawa interactions. In addition to the lower bounds on the mass (m H± ) from CERN 



LEP and Tevatron direct searches 



16 



171 ]. m H ± can also be constrained by the 5-meson 



1 The type-X (type-YjTHDM is referred to as type-IV (type-Ill) THDM in Ref. and the type-I' (type- 



II') THDM in Ref. 12). Sometimes the most general THDM, in which each fermion couples to both Higgs 



doublet fields, is called the type-Ill 14[. To avoid confusion we just call them as these type-X and type-Y 
THDMs. 



3 



decay data such as B — > X s ^ 



19 



20 



2l| and B — > rv 22j, |23| , depending on the model 



of Yukawa interaction. The B — > X s 7 results give a severe lower bound, rn H ± > 295 GeV, at 



the next-to-next-to- 
the type-Y THDM 



eading order (NNLO) in the (non-supersymmetric) type-II THDM and 
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21| , but provide no effective bound in the type-I (type-X) THDM 



for tan/3 > 2, where tan/5 is the ratio of the vacuum expectation values (VEVs) of the 
CP-even Higgs bosons. We also discuss the experimental bounds on the charged Higgs 
sector from pur ely leptonic observables r — > fiVi> 24j and the muon anomalous magnetic 



moment 



25 



26). 



We finally discuss the possibility of discriminating between the types of Yukawa inter- 
actions at the LHC and also at the ILC We mainly study collider phenomenology in the 
type-X THDM in the light extra Higgs boson scenario, and see differences from the results in 
the MSSM (the type-II THDM). We discuss the signal of neutral and charged Higgs bosons 
at the LHC, which may be useful to distinguish the type of the Yukawa interaction. The 
feasibility of the direct production processes from gluon fusion gg — > A (H) and the associ- 
ated production from pp — > bbA (bbH) is studied and the difference in the signal significance 
of their leptonic decay channels is evaluated in the type-X THDM and the MSSM. We also 
consider the Higgs boson pair production pp — > AH ± , HH ± , AH and find that the leptonic 
decay modes are also useful to explore the type of the Yukawa interaction. At the ILC, the 
process of e + e~ — > AH is useful to examine the type-X THDM, because the final states are 
completely different from the case of the MSSM. 

In Sec. II, we give a brief review of the types of the Yukawa interactions in the THDM. In 
Sec. Ill, the decay widths and the branching ratios are evaluated in the four different types 
of the Yukawa interactions. Section IV is devoted to a discussion of current experimental 
constraints on the THDM in each type of the Yukawa interaction. In Sec. V, the possibility 
of discriminating the type of the Yukawa interaction at the LHC and the ILC is discussed. 
Conclusions are given in Sec. VI. The formulae of the decay rates of the Higgs bosons are 
listed in the Appendix. 

II. TWO HIGGS DOUBLET MODELS UNDER THE Z 2 SYMMETRY 



In the THDM with isospin doublet scalar fields $i and $2 with a hypercharge of Y = 1/2, 
the discrete Z 2 symmetry (<£>! — > $x and $2 "~ > — ^2) m ay be imposed to avoid FCNC at 
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TABLE I: Variation in charge assignments of the Z 2 symmetry. 



the lowest order [lCj . The most general Yukawa interaction under the Z 2 symmetry can be 
written as 



/-THDM 

•^"yukawa 



Q L YJ& u u R - Q L Y d $ d d R - L L Y e $ e £ R + H.c, 



(1) 



where $j (/ = u, d or £) is either $x or $ 2 - There are four inde pen dent Z 2 charge assignments 
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12j. In the type-I THDM, 



on quarks and charged leptons, as summarized in TABLE [T] 
all quarks and charged leptons obtain their masses from the VEV of $2- In the type-II 
THDM, masses of up- type quarks are generated by the VEV of $2> while those of down- 
type quarks and charged leptons are acquired by that of $ x . The Higgs sector of the MSSM 
is a special THDM whose Yukawa interaction is of type-II. The type-X Yukawa interaction 
(all quarks couple to $ 2 while charged leptons couple to $1) is used in the model in Ref. 
The remaining one is referred to as the type-Y THDM. 

The most general Higgs potential under the softly broken Z 2 symmetry is given by 

Ai , _ + T \ o As 



yTHDM = m 2 $ t $i + m 2 $ t $2 _ m 2 + + ^t^a + ^($t$ 2 )^ 



+ A 3 ($l$i)($ 2 $ 2 ) + A 4 ($ T 1 $ 2 )(^$ 1 ) + 



it, 



A, 



2 L 



($1$ 2 ) 2 + (^$i) 2 



(2) 



where the parameters ml and A5 are complex, in general. In this paper, they are taken to 
be real by assuming CP invariance. The Higgs doublet fields can be parametrized as 



■j=(vi + hi -izi 



and the mass eigenstates are defined by 







(3) 




R(/3) 




(4) 
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TABLE II: The mixing factors in Yukawa interactions in Eq. © 
where the rotation matrix is given by 

(cos 9 — sin 9 \ 
(5) 
sin 9 cos 9 I 

There are five physical Higgs bosons, i.e., two CP-even states h and H, one CP-odd state 
A, and a pair of charged states H ± , and z and u^ 1 are Nambu-Goldstone bosons that are 
eaten as the longitudinal components of the massive gauge bosons. The eight parameters 
mf-mj and Ai— A5 in the Higgs sector are replaced by eight physical parameters: i.e., the 
VEV v = y/v \ + v 2 — 246 GeV, the mixing angles a and P (tan/3 = v 2 /vi), the Higgs boson 



masses m h , m H , m A , m H ±, and the soft breaking mass parameter M = m 3 /y/sm P cos p. The 
mixing angle a is defined such that h is the SM-like Higgs boson when sin(/3 — a) = 1. 
The Yukawa interactions are expressed in terms of mass eigenstates of the Higgs bosons 

as 



/-•THDM 

•^"yukawa 



f=u,d,i 

\ V2V ud _ 
I V 



(m u C A P L + mt&Pn) dH + + ^*& u 



L*-R 



H + + H.c. 



(6) 



where Pl/r are projection operators for left- /right-handed fermions, and the factors ^ are 
listed in TABLE [Hj 

For the successful electroweak symmetry breaking, a combination of quartic coupling 
constants should satisfy the condition of vacuum stability 
bounds from perturbative unitarity 



21, 



28(. We also take into account 



to restrict parameters in the Higgs potential [301.131]. 



The top and bottom Yukawa coupling constants cannot be taken too large. The requirement 
\Y t ^\ 2 < 7r at the tree level can provide a milder constraint 0.4 < tan/3 < 91, where 
\Y t \ = (y/2 /v)m t cot P and \Y b \ = (V2/v)m b tan/5. 
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FIG. 1: Total decay widths of H, A and i7 ± in the four different types of THDM as a function 



of the decaying scalar boson mass for several values of tan/3 under the assumption 



m , 



m H± and M 



1 GeV. The SM-like limit sin(/3 — a) = 1 is taken, where h is the 



SM-like Higgs boson. 

III. DECAYS OF HIGGS BOSONS 



In this section, we discuss the difference in decays of the Higgs bosons for the types of 
Yukawa interactions in the THDM. We calculate the decay rates of the Higgs bosons and 
evaluate the total widths and the branching ratios. In particular, we show the result with 



sin(/5 — a) 



32 



331 ] . where h is the SM-like Higgs boson while the VEV of H is very 
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FIG. 2: Decay branching ratios of H, A and in the four different types of THDM as a function 
of tan/3 for m H = m A = m H± = 150 GeV and M = 149 GeV. The SM-like limit sin(/3 — a) = 1 is 
taken, where h is the SM-like Higgs boson. 



small or negligible. The decay pattern of h is almost the same as that of the SM Higgs 
boson with the same mass at the leading order except for the loop-induced channels when 
sin(/3 — a) = 1. In this case, H does not decay into the gauge boson pair at tree level, so 
it mainly decays into fermion pairs 2 . We note that A and H ± do not decay into the gauge 
boson pair at the tree level for all values of sin(/5 — a). 

The decay patterns are therefore completely different among the different types of Yukawa 



interactions 
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121 ]. For the decays of H and A, we take into account the decay channels 



2 In the case with a more complicated mass spectrum a heavy Higgs boson can decay into the states which 
contain lighter Higgs bosons [34 1 . 
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FIG. 3: Decay branching ratios of h and H in the four different types of THDM as a function of 
tan (3 for m h = 120 GeV, m H = 150 GeV, M = 148 GeV and sin 2 (/3 - a) = 0.96. 



of qq, £ + £~ , (WVK^, ZZ^*>) at the tree level, and gg, 77, Z7 at the leading order, where q 
represents s, c, b (and t), and I represents fi and r. Running masses for b, c, and s quarks 
are fixed as nib = 3.0 GeV, WT C = 0.81 GeV and mj = 0.046 GeV, respectively. For the decay 
of the charged Higgs boson, the modes into tb, cb, cs, tu, and fiu are taken into account as 
long as they are kinematically allowed. The analytic formulas of each decay rate are given 
in the Appendix for completeness. They are consistent with the previous results for the 
type-II THDM (2]. 

In FIG. [U the total widths of H , A and H are shown as a function of the mass of decaying 
Higgs bosons for several values of tan f3 in the four different types of Yukawa interactions. 
We assume sin(/3 — a) = 1 and m§ = m H = m A = m H ± . To keep perturbative unitarity for 



a wide region of tan/3, the soft breaking parameter is taken to be M = m $ — 1 GeV 30]. 
The widths strongly depend on the types of Yukawa interactions for each tan j3 value before 
and after the threshold of the ti (tb) decay mode opens. 

In FIG. [U the decay branching ratios of H, A and H ± are shown for m $ = 150 GeV, 
sin(/3 — a) = 1, and M = m$ — 1 GeV as a function of tan/3. In the type-I THDM, the 
decay of H into a gauge boson pair 77 or can increase for large tan f3 values, because all 



9 



the other fermionic decays (including the gg mode) are suppressed but the charged scalar 
loop contribution to 77 and Z'j decay modes is not always suppressed for large tan (3. Such 
an enhancement of the bosonic decay modes cannot be seen in the decay of A since there 
is no AH + H~ coupling. In the type-X THDM, the main decay mode of H and A is t + t~ 
for tan/3 > 2, and the leptonic decays t + t~ and become almost 99% and 0.35% for 

tan f3 > 10, while the bb (or gg) mode is always the main decay mode in all other cases. In 
the type-Y THDM, the leptonic decay modes of H and A are rapidly suppressed for large 
tan/? values, and only the branching ratios of bb and gg modes are sizable for tan/? > 10. 
In charged Higgs boson decays with m H± = 150 GeV, the decay into rv is dominant in the 
type-I THDM, the type-II THDM and the type-X THDMs for tan/3 > 1, while hadronic 
decay modes can also be dominant in the type-Y THDM jll| . 

In FIG. [3j we show the decay branching ratios of CP-even Higgs bosons, where = 120 
GeV, m<j,(= m H = m A = m H± ) = 150 GeV, M = 148 GeV, and sin 2 (/3 - a) = 0.96 are 
taken. Because of the CP-even Higgs boson mixing, the lightest Higgs boson h is no longer 
purely SM-like. Instead, H can decay into massive gauge bosons via off-shell modes such as 
WW* and ZZ*. Decay patterns of h and H depend on tan /3 and also on the type of Yukawa 
interaction. When tan/3 ~ 5, the angle a is nearly zero. In such case coupling constants 
become small, so that some of fermionic decay modes are suppressed. In order to satisfy 
the unitarity constraints for the large tan f3 region, the soft breaking mass scale M must be 



degenerate to the mass of the decaying bosons [30 ] . 



IV. CONSTRAINTS FROM THE CURRENT EXPERIMENTAL DATA ON 
THDMS 

One of the direct signal of the THDM is the discovery of extra Higgs bosons, which have 



been searched at the LEP and Tevatron 
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171 ]. Indirect contributions of Higgs bosons to 



precisely measurable observables can also be used to constrain Higgs boson parameters. In 
this section, we summarize these experimental bounds. 

A direct mass bound is given from the LEP direct search results CIS TTijjQ > 92.8 GeV 
for CP-even Higgs bosons and m A > 93.4 GeV for CP-odd Higgs bosons in supersymmetric 
models. The bound for charged Higgs boson has also been set as m H ± > 79.3 GeV [161 ] . 

In the THDM, one-loop contributions of scalar loop diagrams to the rho parameter are 
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expressed as 

<Vthdm = ^ Gf 7^2 { F A{m 2 A ,m 2 H ±) 

-c 2 a _p [F A {m 2 h ,m 2 A ) - F A {m 2 h ,m 2 H± )] - s 2 a _ p [F A {m 2 H , m 2 A ) - F A (m 2 H , m 2 H± )] } , (7) 

where F A (x, y) = ^(x+y) — In | with F A (x, x) = 3 . These quadratic mass contributions 
can cancel out when Higgs boson masses satisfy the following relation: (i) m A ~ m H± , (ii) 
m H ~ m H± with sin(/3 — a) ~ 1, and (iii) m h ~ m„ + with sin(/5 — a) ~ 0. These relations 



correspond to the custodial symmetry invariance 35|, [36j . This constraint is independent of 
the type of Yukawa interaction. 

When m H± < m t — m^, the top quark can decay into the charged Higgs boson. The decay 
mode t — > H + b has been studied by using the Tevatron data [37J]. The partial decay width 
is calculated by using the factors in TABLE [Til as 



x < m t 



+ 4m?m^d !> . (8) 



where X(x, y) is defined in Eq. flA5p . The higher order corrections can be found in Ref. 38]. 
The decay branching fraction strongly depends on the type of Yukawa interaction. In the 
type-I (type-X) THDM, the decay mode can be sizable only for small tan /3. On the other 
hand, it can also be substantial for the large tan/5 region in the type-II (type-Y) THDM 
since the bottom quark Yukawa coupling receives considerable enhancement. This fact gives 
a lower bound 1 < tan (3 for THDMs and an upper bound tan (3 < 60 for the type-II (type-Y) 
THDM below m H± < 130GeV [37]. 

It has been known that the charged Higgs boson mass in the type-II THDM is stringently 
constrained by the precision measurements of the radiative decay of b — > 57 at Belle [39| and 
BABAR |40j as well as CLEO [411 ] . The process b — >• 37 receives contributions from the W 
boson loop and the charged Higgs boson loop in the THDM. A notable point is that these 
two contributions always work constructively in the type-II (type-Y) THDM, while this is 



3 There are other (relatively smaller in most of the parameter space) contributions to the rho parameter 
in the THDM, i.e., those from the diagrams where the gauge boson (as well as Nambu-Goldstone boson) 
and the Higgs boson are running together in the loop[2[. We have included these effect in our numerical 
analysis. 
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FIG. 4: Predictions of the decay branching ratio for b — > 57 are shown at the NLO approximation 
as a function of rn H± and tan (3. The dark (light) shaded band represents la (2a) allowed region 
of current experimental data. In the left panel, solid (dashed) curves denote the prediction for 
tan/3 = 2 (50) in various THDMs. In the right panel, solid, dashed and dot-dashed curves are 
those for m H ± = 100, 300 and 600 GeV, respectively. 



not the case in the type-I (type-X) THDM llj. In FIG. HI we show the branching ratio 
of B — > Xg'-f for each type of THDM as a function of m H± (left-panel) and tan/? (right- 
panel), which are evaluated at the next-to-leading order (NLO) following the formulas in 
Ref. 181 ] . The SM prediction at the NLO is also shown for comparison. The theoretical 
uncertainty is about 15% 4 in the branching ratio ( as indicated by dotted curves in FIG. H), 
which mainly comes from the pole mass of charm quark m[! ole = 1.65 ± 0.18 GeV [42]. 
The experimental bounds of the branching ratio are also indicated, where the current world 
average value is given by B{B -> X s7 ) = (3.52 ±0.23 ±0.09) x 10" 4 [42]. It is seen in FIG. HI 
that the branching ratio in the type-I (type-X) THDM lies within the 2 a experimental 
error in all the regions of m#± indicated for tan f3 > 2, while that in the type-II (type-Y) 
THDM is far from the value indicated by the data for a light charged Higgs boson region 
(m H± < 200 GeV). In the right figure, a cancellation occurs in the type-I (type-X) THDM 
since there are destructive interferences between the W boson and the H ± contributions. It 



In Ref. 



18 1, the theoretical uncertainty is smaller because the value for the error in m^ ole /ml° le is taken 



to be 7%, which gives main uncertainty in the branching ratio. 
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FIG. 5: Bounds from B —* ri/ (left panel) and tau leptonic decay (right panel) on m-^± as a 
function of tan/3 are shown. The dark (light) shaded region corresponds to the 2a (la) exclusion 
of these experimental results. In the type-II THDM the wide parameter space is constrained by 
B — > rv, while only the tau leptonic decays are important for the type-X THDM. 

is emphasized that the charged Higgs boson could be light in the type-I (type-X) THDM 
under the constraint from B — ► Xgj results. We note that in the MSSM the chargino 
contribution can compensate the charged Higgs boson contribution |44j. This cancellation 
weakens the limit on m H± from b — > 57 in the type-II THDM, and allows a light charged 
Higgs boson as in the type-I (type-X) THDM. 

We give some comments on the NNLO analysis, although it is basically out of the scope 
of this paper. At the NNLO, the branching ratio for b —>■ sj has been evaluated in the 



SM in Ref. 



2Cj, 



21] . The predicted value at the NNLO is less than that at the NLO 



approximation in a wide range of renormalization scale. In Ref. [20j, the SM branching 
ratio is (3. 15 ±0.23) x 10~ 4 , and the lower bound of the m H± , after adding the NLO charged 
Higgs contribution, is estimated > 295 GeV (95% CL) in the type-II (type-Y) 

THDM [20] 5 . On the other hand, in the type-I (type-X) THDM, although the branching 
ratio becomes smaller as compared to the NLO evaluation, no serious bound on m H± can 
be found for tan/5 > 2. Therefore, charged Higgs boson mass is not expected to be strongly 



5 In Ref. [21[ the NNLO branching ratio in the SM is calculated as (2.98 ±0.26) x 10" 4 , and the mass bound 
is a little bit relaxed. 
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constrained in the type-I (type-X) THDM even at the NNLO, and our main conclusion that 
the type-I (type-X) THDM is favoured for m H± < 200 GeV based on the NLO analysis 
should not be changed. 

The decay B — > tv has been discussed in the type-II THDM 
B{B + -> t + v t ) = (1.4 ± 0.4) x 10" 4 are obtained at the B factories 
branching ratio can be written as 22j, |46| 



22 



231 ] . The data for 
45]. The decay 



42 



B(B+ 



T IV J THDM 



ni 



B 



m 



pdpt 

SAW 



(9) 



In FIG. 15] the allowed region from the B — > tv results is shown in the type-II THDM. The 
dark (light) shaded region denotes the 2a (la) exclusion from the B — > tv measurements. 
The process is important only in the type-II THDM with large tan (3 values. The other types 
of Yukawa interactions do not receive constraints form this process. 

The rate for the leptonic decay of the tau lepton r — > uvv can be deviated from the 
SM value by the presence of a light charged Higgs boson [24|. The partial decay rate is 
approximately expressed as [221 ] 



-pTHDM 
■pSM 



2m 



m 



+ 



4m 4 H± 



A 5 



(10) 



where k(x) = g{x)/f(x) is defined by f(x) = 1 — 8x — 12x 2 In x + 8x 3 — x A , and g(x) = 
1 + 9a; — 9a; 2 — x 3 + 6a; (1 + x) In a;. In the type-II (type-X) THDM, the leptonic Yukawa 
interaction can be enhanced in the large tan/? region. Hence, both the models are weakly 
constrained by tau decay data, as in FIG. [5] 

The precision measurement for the muon anomalous magnetic moment can give mass 
bound on the Higgs boson in the SM [25j. This constraint can be applied for more general 
interaction, including THDMs [26]. At the one-loop level, the contribution is given by 



4tt 2 V2 



E 



tl 2 



In 



i>°=h,H "V 



m. 



m 



In 



777,4 11 

— H 

777 2 , 6 



6m 2 H± 



This process is also purely leptonic and only gives milder bounds on the Higgs boson masses 
for very large tan/? values in the type-II (type-X) THDM. It gives no effective bound on 
the type-I (type-Y) THDM. It is also known that the two-loop (Barr-Zee type) diagram can 

. The contribution can be large because of the enhancement 



significantly affect a. 



47. 
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factors of m|/m^ and also of the mixing factors ^ as 



BZ NfQ 2 f G F am 2 b 



47T 3 f 2 



where 



/(*) 



1 — 1x{\ — x) x(l — x) 
In 



— / G?X " 

2 J x(l-x)-z 

1 1 
c&e^ : In 



z 

x(l — x) 



(12) 



(13) 
(14) 



2 J x(1 — x) — z z 
The contribution from this kind of diagram is only important for large tan/? values with 
smaller Higgs boson masses in the type-II THDM. For the other types of THDM, it would 
give a much less effective bound on the parameter space. 



V. COLLIDER SIGNALS IN THE TYPE-X THDM AT THE LHC AND THE ILC 



We discuss the collider phenomenology of the models at the LHC and the ILC. There 
have already been many studies on the production and decays of the Higgs bosons in the 
type-II THDM, especially in the context of the MSSM, while the phenomenology of the other 
types of THDMs has not yet been studied sufficiently Recently, the type-X THDM has been 
introduced in the model to explain phenomena such as neutrino masses, dark matter, and 
baryogenesis at the TeV scale We therefore concentrate on the collider signals in the 
type-X THDM, and discuss how we can distinguish the model from the type-II THDM (the 
MSSM), mainly in scenarios with a light charged Higgs boson (100 GeV < m H ± ~ 300 
GeV). (Such a light charged Higgs boson is severely constrained by the b —>■ sj result in 
the non-supersymmetric type-II THDM and the type-Y THDM, while it can be allowed in 
the MSSM and the type-X (type-I) THDM.) As we are interested in the differences in the 
types of the Yukawa interactions, we focus here on the case of m H± ~ m A ~ m H with 
sin(/5 — a) ~ 1 for definiteness. 



A. Charged Higgs boson searches at the LHC 

A light charged Higgs boson with m H± < mt — mt can be produced in the decay of 
top quarks at the LHC. The discovery potential for the charged Higgs boson via the ti 
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production has been studied in the MSSM [37[. Assuming an integrated luminosity of 30 
fb -1 , the expected signal significance of the event tt — > H ± W Zf bb — > tvTv T bb is greater 
than 5<T for m H ± < 130 GeV for tan f3 < 2 and tan (3 > 20 37|] . The same analysis can 
also be applied for the type-X THDM, in which a similar number of H can be produced 
when tan/5 ~ 0(1)- The main decay mode (tu) is common in the type-II THDM and the 
type-X THDM, except for very low tan (3 values. Thus, searching for neutral Higgs bosons 
is important to distinguish the difference in the types of Yukawa interactions. 



For m H ± > m t , charged Higgs bosons can be produced in qq/gg — > tbH , gb — > tH 



50|, gg (qq) - H+H~ 



51 



49, 



52j | and gg (bb) — > H ± W T [53J. These processes, except for 



the H + H~ production, are via the Yukawa coupling of tbH", so that the cross sections 
are significant for tan/3 ~ 0(1) or tan/5 > 10-20 in the type-II THDM and only for 
tan/3 ~ 0(1) in the type-X THDM. The type of Yukawa interaction in the THDM can 
then be discriminated by measuring the difference in decay branching ratios of H ± . In the 
type-II THDM H ± mainly decay into tb, while rv is dominant for tan (3 > 10 in the type-X 
THDM. 



B. Neutral Higgs boson (A and H) production at the LHC 



At the LHC, the type of the Yukawa interaction may be determined in the search for 



neutral Higgs bosons through the direct production via gluon fusion gg 
vector boson 



bbA (bbH) 



57|, 



A/H |52|, |54j, 



usion V*V* — > H, V = W, Z [55L l56j]. and also via associated production pp 



50 



58 



59J can also be useful for 



The production process pp —> ttA (tiH) 
tan (3 ~ 1. We discuss the possibility of discriminating the type of the Yukawa interaction by 
using the production and decay processes of A and H for sin(/3 — a) = 1. Additional neutral 
Higgs bosons A and H are directly produced by the gluon fusion mechanism at the one-loop 
level. When sin(/? — a) ~ 1, the production rate can be significant due to the top quark 
loop contribution for tan/? ~ 1 and, in the MSSM (the type-II THDM), also for large tan/3 
via the bottom quark loop contributions fo^l ]. Notice that there is no rate of V*^* — > A 
because there is no VVA coupling, and that the production of H from the vector boson 
fusion V*V* — > H is relatively unimportant when sin(/3 — a) ~ 1. The associate production 
process pp — » bbA (bbH) can be significant for large tan (3 values in the MSSM (the type-II 
THDM) [5J 
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In the MSSM (the type-II THDM), the produced A and H in these processes decay 
mainly into bb when sin(/3 — a) ~ 1, which would be challenging to detect because of huge 
QCD backgrounds. Instead, the decays into a lepton pair t + t~ would be useful 

for searches of A (and H). However, the decay branching ratios of A — > t + t~ pr) are 
less than 0.1 (0.0004). A simulation study 60] shows that the Higgs boson search via the 
associate production bbA (bbH) is better than that via the direct production from gluon 
fusion to see both r + r~ and [i + pT modes, especially in the large tan/3 area. The largest 
background is the Drell-Yan process from 7* jZ* — f t + t~ (and p + p~). The other ones, such 
as tt, bb and W + jet, also contribute to the backgrounds. The rate of the r + r~ process 
from the signal is much larger than that of the p + p~ one. However, the resolution for tau 
leptons is much broader than that for muons, so that for relatively small itla {ttih) the p + p~ 
mode can be more useful than the r + r" mode 61]. 

In the type-X THDM, signals from the associate production pp — > bbA are very difficult to 
detect. The production cross section is at most 150 fb for m A = 150 GeV at tan/3 =1 60], 
where the branching ratio A — > t + t~ and A — > yU + /i~ are small, and the produced signals 
are less for larger values of tan (3. On the other hand, the direct production from gg — > Aj H 
can be used to see the signal. The cross sections are significant for tan/5 ~ 1, and they 
decrease for larger values of tan/3 by a factor of 1/ tan 2 /3. However, the branching ratios of 
A/H — > t + t~ dominate over those of A/H — > bb for tan/3 > 2 and become almost 100% for 
tan/5 > 4 (see FIG. |2]). Therefore, large significances can be expected around tan/5 ~ 2 in 
the type-X THDM. 

In FIG. El we show the expected signal significance of the direct production pp(gg) — > 
A/H — ► t + t~ in the type-X THDM and the MSSM at the LHC, assuming an integrated 
luminosity of 30 fb -1 . The mass of the CP-odd Higgs boson m A is taken to be 150 GeV 
in both models, while m H is 150 GeV for the results of the type-X THDM and that in the 



MSSM is deduced using the MSSM mass relation. For the detai 



ed analysis of background 



simulation, we employed the one shown in the ATLAS TDR 60]. The basic cuts, such as 
the high p T cut and the standard A/H — > r + r~ reconstruction, are assumed 6^]. We can 
see that, for the search of the direct production, the signal significance in the type-X THDM 
can be larger than that in the MSSM for tan/5 < 5. In particular, the signal in the type-X 
THDM can be expected to be detectable (S/yB > 5) when tan/5 ~ 2 for the luminosity 
of 30 ftr 1 . For smaller values of uia (jnu), the production cross section becomes large so 
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FIG. 6: Signal significance (S/VB) for gluon fusion gg — > A/H (left panel) and the associated 
production pp — > bbA (bbH) (right panel) with the t + t~ final state in the type-X THDM and the 
MSSM. In both figures, the dashed and the solid curves represents the expected values of signal 
significance for A/H and summed over A and H. The red (thin) curves denote the significance in 
the type-X THDM, while the black (thick) solid curves denote that in the MSSM. 



that the signal rate is more significant, but the separation from the Drell-Yan background 
becomes more difficult because the resolution of the tau lepton is broad. Therefore, the 
significance becomes worse for m A (m H ) < 130 GeV. 

When A and H are lighter than 130 GeV, the [i + mode can be more useful than 
the t + t~ mode. The resolution of muons is much better than that of tau leptons, so that 
the invariant mass cut is very effective in reducing the background from 7*/Z* — > /i + pT . 

A/H — > has been studied in the SM and the 

, |6j|. We evaluate the si gnal significances of gg — > A/H — > fi + [i~ in the 



The feasibility of t he p rocess gg 



MSSM in Refs. 



60 



type-X THDM by using the result in Ref . 6j| . In TABLE IHIl we list the results of the 
significance in the SM and the type-X THDM. According to Tao Han's paper, the basic 



kinematic cuts of p T > 20 GeV, \rj\ < 2.5 and the invariant mass cuts as m 



A/H 



2.24 GeV 



< M < rn A/H + 2.24 GeV are used 6 . The integrated luminosity is assumed to be 300 



This choice for the invariant mass cut is rather severe; i.e., it requires the precise determination of m A and 
ni H . If the range of the invariant mass cut is taken to be double, roughly speaking, background events 
also become double. This would suppress the signal significance in TABLE IlIII by a factor of ~ 1/V2. 
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m 6 [GeVl 


SM (H S m) 


MSSM(A) 


Type-X (H) 


Type-X (A) 


Type-XfSum) (m A ~ rriij) 


115 


2.41 


1.31 


4.31 


12.0 


16.3 


120 


2.51 


1.49 


4.89 


13.4 


18.3 


130 


2.25 


1.81 


5.78 


15.6 


21.4 


140 


1.61 


2.11 


6.60 


17.5 


24.1 



TABLE III: Expected signal significances for gg — > (f> — > (0 = H^m, H, A) in the SM, the 

MSSM, and the type-X THDM. For the results in the MSSM tan (3 = 2 is taken, and for that in 
the type-X THDM tan (3 = 2 and sin(/3 — a) = 1 are taken. The integrated luminosity is assumed 
to be 300 ft) -1 . 

fb _1 . For the results in the type-X THDM, tan/? = 2 and sin(/3 — a) = 1 are taken. The 
results show that the significance can be substantial for m A > 115 GeV when tan (3 = 2. 
For smaller masses of the extra Higgs bosons, the cross section for the signal processes can 
be larger, but the invariant mass cut becomes less effective in the reduction of the Drell-Yan 
background because of the smaller mass difference between m z and m A ; hence, the signal 
significance becomes worse. The tan/? dependence in the signal significance for the muon 
final states is also shown in FIG. [71 The shape of the curves is similar to that for the tau 
lepton final state in FIG [HI 

In summary, we would be able to distinguish the type-X THDM from the MSSM by 
measuring the leptonic decays of the additional Higgs bosons produced via the direct search 
processes gg — > A/ H — » t + t~ and the associated processes pp — > bbA (bbH). First, 

if a light scalar boson is found via gg h — > 77 or W + W~ — > h — > 77 (t + t - ) and if 
the event number is consistent with the prediction in the SM, then we know that the scalar 
boson is of the SM or at least SM-like: in the THDM framework this means sin(/3 — a) ~ 1, 
assuming that it is the lightest one. Second, under the situation above, when the associated 
production pp — ■> bbA (bbH) is detected at a different invariant mass than the mass of the 
SM-like one and no gg — > A/H — > t + t~ (fi + pr) is found at that point, we would be able 
to identify the MSSM Higgs sector (or the type-II THDM) with high tan /3 values. On the 
other hand, the type-X THDM with low tan (3 would be identified by finding the signal from 
the gluon fusion process without that from the bbr + r~ . If signals from the direct production 
processes are found but the number is not sufficient, then the value of tan/5 would be around 
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gg^ 




FIG. 7: Signal significance (S/^/B) for the final state from gluon fusion gg — > A in the 

type-X THDM and the MSSM. The red (thin) curves denote significance in the type-X THDM, 
while the black (thick) solid curves denote that in the MSSM. 

6-10 (mt cot (3 ~ m\, tan ft), where the rates in the MSSM and the type-X THDM are similar. 
In this case, it would be difficult to distinguish these models from the above processes. As 
we discussed in the next subsection, Higgs pair production processes pp — > AH ± , HH ± , and 
AH can be useful to measure the Yukawa interaction through branching fractions, because 
these production mechanism do not depend on tan (3 in such a situation. 

We have concentrated on sin(/3 — a) ~ 1 in this analysis because the parameter is moti- 
vated in Ref. 7|]. Here we comment on the situation without the condition sin(/3 — a) = 1. 
If sin(/5 — a) is not close to unity, our conclusion can be modified. The production cross sec- 
tions of gg — > A/H and pp — > bbA(bbH) can be enhanced in the type-X THDM for tan (3 > 1 
since the factor (sin a/ sin (3) of quark-Higgs couplings can be large in a specific region of 
the parameter space. The signal of the CP odd Higgs boson A can then be significant. On 
the other hand, the CP even Higgs boson H can decay into WW* when sin(/5 — a) is not 
unity. This would reduce leptonic branching fractions. The signal can be enhanced only for 
large tan (3 regions because leptonic decays are significant only for such a parameter space. 
We also note that H can be produced significantly by the vector boson fusion mechanism 
in a mixing case. 
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FIG. 8: Production cross sections of pp — > AH^ and AH are shown at the leading order as a 
function of scalar boson masses in the THDM where m$ = m H = m A = rn H± are chosen. The long- 
dashed, dashed and solid curves denote the pair production of AH + , AH~ and AH, respectively. 
The cross section of pp — ► HH^ is the same as those of pp — > AH^ when sin(/3 — a) = 1. 

C. Pair production of extra Higgs bosons at the LHC 

The types of the Yukawa interactions can be studied using the process of qq' — > — > 
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63 



64 



651 ] and gg — > Z* — > Aff [52j |. unless the extra Higgs bosons 



AH* {HH ± ) (52 

if, A and/or if* are too heavy 7 . Hadronic cross sections for these processes are shown at 
the leading order in FIG. [8] as a function of the mass of the produced scalar boson $, where 
m $ = m H = m A = m H ±. Expected rates of AH ± (sum) and AH are about 143 fb and 85 
fb for m H± = 130 GeV and about 85 fb and 50 fb for m H± = 150 GeV, respectively. The 
NLO QCD corrections are expected to enhance these rates typically by about 20% [62j, |64 |. 
The production rates are common in all the types of THDMs, because the cross sections 
are determined only by the masses of the produced scalar bosons. Therefore, they are very 
sensitive to the difference in the decay branching ratios of the produced Higgs bosons. In 
the MSSM, the bbr v (bhr + T~) events can be the main signal of the AH ± and HH ± {AH) 
processes, while in the type-X THDM (tan/3 > 2), t+t^t^v (t + t~t + t~) would be the main 



7 When the mixing between h and H is large, the hH ± production can also be important 65 1. 
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signal events. 

In the MSSM, the parton level b ackg round analysis for the AH (HH^) production 



process has been performed in Ref. 



64 



65| by using the bbrv final state. The largest 



background comes from qq' — > Wg* — > Wbb, which can be reduced by basic kinematic cuts 
and the invariant mass cut of bb, as well as by the kinematic cut to select hard hadrons 
from the parent r ± from H ± . It has been shown that a sufficient signal significance can be 



obtained for smaller masses of Higgs bosons 64 ] . 



In the type-X THDM, the produced AH* {HH*) and AH pairs can be studied via the 
leptonic decays. Hence these channels can be useful to discriminate the type-X THDM from 
the MSSM. Assuming an integrated luminosity of 300 fb _1 , 8.6 x 10 4 and 5.1 x 10 4 of the 
signal events are produced from both AH and HH ± production for m $ = 130 GeV and 
150 GeV, respectively, where = m H = m A = m H± . A and H (H ) decay into t + t~ 
(tv) by more than 95% and 95% (99%) for tan/5 > 4, respectively. The purely leptonic 
signal would have an advantage in the signal to background ratio because the background 
from the intermediate state qq' — > Wg* would be negligible. For tan/5 = 7, the produced 
AH and HH ± pairs almost all (more than 99%) go to t + t~tv final states. The numbers 
of signal and background are summarized in TABLE IIV1 The signal to background ratio for 
t + t~tv final state is not so small 0(0.1-1), before cuts 8 . The backgrounds are expected 
to be reduced by using high-p T cuts, hard hadrons from the parent tau leptons from H ± , 
and invariant mass cuts for t + t~ from A and H, in addition to the basic cuts. However, 
the signal significance strongly depend on the rate of miss identification of hadrons as tau 
leptons, so that a realistic simulation is necessary. 

We also would be able to use the fi + ^Tt + v events to identify the AH + and HH + produc- 
tion in the type-X THDM, by using the much better resolution of in performing the 
invariant mass cuts. For 300 fb _1 , the AH + and HH + process can produce about 0(100) 
of [i + fi~T + b> events for m A = m H = 130 GeV. The number of background events is about 
3.1 x 10 5 of fi + fi~T + v from ZW ± production. Signals and background for yU + /i _ r + z/ events 
are also summarized in TABLE IIVI The background can be expected to be reduced by 
imposing a selection of the events around the invariant mass of m A ~ and the high p T 



The production may give a much larger cross section for background events. It may also be reduced 

considerably by kinematic cuts. 
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AH ± ,HH ± {m^ = 130 GeV) 


AH ± ,HH ± {m^ = 150 GeV) 


ZW ± 


T + T TV 


8.4 x 10 4 


5.0 x 10 4 


3.2 x 10 4 


/1 + [X~TV 


3.0 x 10 2 


1.8 x 10 2 


3.1 x 10 4 



TABLE IV: Events for the t + t~tv and [i + [i~tv final states from the Higgs boson pair production 
and ZW^ background. The signal events are summed over AH^ and HH^ 1 . The integrated 
luminosity is taken to be 300 fb _1 . Values for the decay branching ratios are taken to be B(A/H — > 
t + t~) = 0.99, B(A/H -> fi + fi-) = 0.0035, and B(H ± -> tv) = 0.99, which correspond to 
the values for tan/3 > 7. The cross section of pp — > ZW^ is evaluated as cr^w = 27.7 pb by 
CalcHEP y. 





Aff(m $ = 130 GeV) 


AH(m^ = 150 GeV) 


zz 




2.5 x 10 5 


1.5 x 10 5 


3.6 x 10 3 




1.8 x 10 2 


1.0 x 10 2 


7.1 x 10 3 



TABLE V: Events for the t + t t + t and r + r fi + fi final states from the Higgs boson pair 
production and ZZ background. The integrated luminosity is taken to be 300 fb _1 . The cross 



section of pp — ► ZZ is evaluated as azz = 10.5 pb by CalcHEP 



66]. 



cuts. Hard hadrons from the decay of r's from H + can also be used to reduce the back- 
ground. In the MSSM, much smaller signals are expected, so that this process can be a 
useful probe of the type-X THDM. 

In a similar way, we may use AH production [52] to identify the type-X THDM. For the 
t + t~t + t~ decay mode, the signal is evaluated approximately as 2.5 x 10 4 events, assuming 
300 ftr 1 for m A = m H = 130 GeV and tan (3 = 7. The main background may come from the 
qq — > ZZ process. We also consider the n + fi~T + r~ decay mode. The number of signal event 
is (9(100) for an integrated luminosity 300 fb _1 . The numbers of signal and background event 
are listed in TABLE [V] It would be valuable to use the detailed background simulation. 
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FIG. 9: The production cross section of e + e 



AH is shown as a function of the Higgs boson 



mass. The dot-dashed, dashed, long-dashed, and solid curves correspond to yfs = 350, 500, 800, 
and 1000 GeV, respectively. 

D. Pair production of extra Higgs bosons at the ILC 

At the ILC, we would be able to test the types of the Yukawa interactions via the pair 



cross section is shown for y/s = 350, 500, 800, and 1000 GeV as a function of m A assuming 
m A = m H in the THDM. The production rate is determined only by m A , and m H at the 
leading order, and is independent of tan/5. (In the MSSM, it depends indirectly on tan/? 
via the mass spectrum.) The signal of the type-X THDM can be identified by searching 
for the events of t + t~t + t~ and h + h~t + t~ . When yfs = 500 GeV, assuming an integrated 
luminosity of 500 fb -1 , the event number is estimated as 1.6 x 10 4 (1.8 x 10 2 ) in the type-X 
(type-II) THDM for t + t~t + t~ , and 1.1 x 10 2 (0.6) for /j, + h~t + t~ assuming m H = m A = 
m H± = 130 GeV, sin(/3 — a) = 1 and tan (3 = 10. This number does not change much for 
tan/3 > 3. The main background comes from the Z pair production, whose rate is about 
4 x 10 2 fb. The event numbers from the background are then 2.3 x 10 2 for t + t~t + t~ and 
4.6 x 10 2 for {i + /j~t + t~. Therefore, the signal should be easily detected in the type-X 
THDM, by which we would be able to distinguish the type-X from the type-II (the MSSM). 

The detailed measurement of the masses of additional Higgs bosons and Yukawa coupling 
constants will make it possible to determine the scenario of physics beyond the SM through 



production of the additional Higgs bosons e + e 



AH [52j, [67[. In Fig. El the production 
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the Higgs physics. 



VI. DISCUSSIONS AND CONCLUSIONS 

We have discussed the phenomenology of the four possible types of Yukawa interactions 
in the THDM with softly broken discrete Z2 symmetry. Although the particle contents are 
the same in these models, their phenomenologies are completely different from each other. In 
some new physics models, two Higgs doublets (plus singlets) are introduced with a definite 
type of Yukawa interaction. Therefore, when additional Higgs bosons are discovered, we may 
be able to distinguish new physics models via the differences between the types of Yukawa 
interactions. 

The differences between the types of the Yukawa interactions largely affect the decay of 
the Higgs bosons. We have evaluated the decay branching ratios in each type of THDM. 
We have then summarized the current experimental bounds on the models from the data of 
b — > S7, B — > tv, r — > jiVu and measurements of muon anomalous magnetic moment. The 
charged Higgs boson mass can be light (m H ± <m t ) in the type-I and type-X THDMs under 
these experimental limits. Such a light charged Higgs boson is also possible in the MSSM. 

We have discussed the collider phenomenology of the THDMs. For definiteness, we have 
concentrated on phenomenological differences between the MSSM Higgs sector and the type- 
X THDM in the case of a light charged Higgs boson. The type-X Yukawa interaction is used 
in TeV-scale models for neutrino masses, dark matter, and baryon asymmetry. At the LHC, 
the type-X THDM can be discriminated from the MSSM by searching for the production 
and decays of the extra Higgs bosons A, H, and H ± , such as gg — > A/H — > £ + £~ , pp — > bbA 
(bbH) — > bb£ + £~ , where £ is r or /i when sin(/5— a) ~ 1. We also discussed the pair production 
processes pp — > AH 1 * 1 , HH ± , and AH to test the type-X THDM. These processes would 
provide distinctive four lepton final states l+i'r^v and £ + £~t + t~ in the type-X THDM, 
while the MSSM Higgs sector can be tested by bbr^v and bbr + T~. Although the realistic 
simulation study is necessary, the type-X THDM can be identified at the LHC with an 
integrated luminosity of 300 fb _1 . At the ILC, the type-X THDM is expected to be studied 
very well by the pair production e + e~ — > AH. The signal should be four leptons (£ + £~t + t~). 

We have discussed phenomenological discrimination of the types of Yukawa interactions 
in the THDM at future experiments. In particular, we have mainly discussed the separation 
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at the LHC and the ILC between the MSSM Higgs sector and the type-X THDM in the 
relatively light charged Higgs boson scenario. By extending this study to include other 
various cases, we can expect that the types of the Yukawa interactions in extended Higgs 
models can be completely determined at current and future collider experiments. Such 
information may be useful to select a true model from many proposed new physics models 
at TeV scales. 
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Note added: 



After this paper was completed, Refs. 
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70j appeared in which similar Yukawa inter- 



action is discussed in a different context. 
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APPENDIX A: HIGGS BOSON DECAY RATES IN THDMS 



In the THDM, the partial decay widths of Higgs bosons decaying into a fermion pair are 
computed at the leading order as 



F(<p -> qq) = N c 
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where the factors £^ are defined in Eq. q = u, d, s, c, t, b; £ — e, /1, r; N c (= 3) is the color 
factor; V u d is the Kobayashi-Maskawa matrix; and 
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Formulas for the decay widths for loop induced decays are given by 
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where fermionic loop functions are given by 

-N c ^f [2 - /3 2 m 2 C (0, 0, m%, m}, m% m 2 )] for <p = h,H 
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with c 



j v ~ 2swc w (^3L ~ ^Qf s vv^-> where T^ L is the third component of the isospin. The 
functions Bq and Co are the Passarino-Veltman functions 7l| . which can be expressed by 
analytic functions as 
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where A denotes the regularized divergence part in the dimensional regularization, and 
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for x < 1 
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The loop functions Iw,Ih±, J\v an d Jh± are common among all types of Yukawa interac- 
tions, which can be found in Ref. 
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